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AU Mic is a young, nearby X-ray active M-dwarf with an edge-on debris disk. Debris disk are the successors of the 
gaseous disks usually surrounding pre-main sequence stars which form after the first few Myrs of their host stars' 
lifetime, when - presumably - also the planet formation takes place. Since X-ray transmission spectroscopy is sensitive 
to the chemical composition of the absorber, features in the stellar spectrum of AU Mic caused by its debris disk can 
in principle be detected. The upper limits we derive from our high resolution Chandra LETGS X-ray spectroscopy are 
on the same order as those from UV absorption measurements, consistent with the idea that AU Mic's debris disk 
possesses an inner hole with only a very low density of sub-micron sized grains or gas. 
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1. Introduction 

The disks around young stars undergo dramatic changes 
during the first ~ 10 Myrs after their host stars' birth, when 
the gas content of the disk largely disappears ([Alexander! 
2008: lMever et al.ll2007tlHernandez et alJl2006lh leaving be- 
hind a so-called debris disk. The Kuiper belt and the aster- 
oid belt are the solar system's analogs of stellar debris disks. 
The main components of an optically thin debris disk are 
small grains with about sub-micrometer sizes, larger bodies 
in the cm range and, possibly, planets, which are thought 
to form in the same time-span. The initial composition of 
the material in the debris disks after the transition phase is 
not well known. Collisions of already formed smaller bod- 
ies replenish the dust in the "older" debris disks, while it is 
not clear, whether this is also the source of the initial dust 
in the debris disk or whether it is remnant proto-planetary 
dust. 



1.1. AU Mic and its activity 

AU Mic is a 12±l Myr old Ml dwarf at a distance of about 
10 pc, which belongs to the /3 Pic moving group (e.g . 
iBarrado v Navascues et all Il999t IZuckerman et alj 12001). 
AU Mic is one of the brightest nearby X-ray emitters 
(\ogLx ~ 29.3) and shows strong flaring activity, making 
AU Mic a valuabl e target for flare stud ies as shown by e.g. 
UV observations (jRobinson et al.ll200ll ). 

At X-ray wavelengths AU Mic has been observed many 
times. The first Chandra observation provided the high- 
est resolution spectrum, but was limited to the wave- 
length range below ~ 25 A (jLinskv et all I2002D . The 
XMM- Newton observation of AU Mic in 2000 was si- 
multaneous with UV and VLA observations rev ealing 
several flares flSmith et al.l 120051: iMitra-Kraev et al.|[2005[ 
iNess et al.ll2003l for line fluxes )■ Furthermore, AU Mic was 
the target of F USE and Hubble Space Telescope STIS ob- 
servations (e.g. iPaeano et al.l [2000; Ro binson et al.l l200lt 



iDel Zanna et al.ll2002f ). aiming at the determination of the 
temperature structure of its chromosphere and corona. 

1.2. AU Mic and its debris disk 

The first indications for cold material around AU Mic go 
back to IRAS data, which exhibit excess emissi on at 60 /im 
(|Mathioudakis fc Dovldll991l;ISong et alj|2002ft . A clear in- 
frared excess at 850 fan in the spectra l energy d i stribu - 
tion (SED) of AU Mic was detected bv ILiu et all (|2004T >. 
clearly pointing to the existence of a debris disk. By assum- 
i ng an optical thin disk at a single temperature, ILiu et aD 
(2004) derived a mass of 0.01 M ffi and a temperature of 
40 K for the disk. These values have been c onfirmed by 



Rebull et al.1 (120081) from Sp itzer data (see also lChen et al.l 
2005h . iMetchev et al.l (|2005l ) also confirmed the dust mass 
of about 0.01 M ffi composed of grains in the submicron 
regime by modelling the opti cal, near-IR and SED data. 

Optical observations bv jKalas et ail |2004), initiated 
shortly after ILiu et all (|2004f) , clearly showed the presence 
of a debris disk around AU Mic with a radial extent of 
at least 210 AU and almost perfectly edge-on. The disk 
has then been subsequently studied at optical wavelengths 
with, e.g., the Hubble Space Telescope (HST) and adap- 
tive optics, making it one of the most well stu died debris 
disks. The models derived bv lKrist et al.1 (|2005|) from their 
HST observation restrict the disk inclination to i > 89°; 
they also c onfirmed th e sma ll-scale brightness variations 
detected bv lKalas et alj (|2004f) . These disk inhomogeneities 
can be readily explained by the existence of orbiting plan- 
ets, how ever, no clear signatures of a planet have been found 
to date (|Hebb et alj|2007t IMetchev et al.l 1 2 0051) . 

Two studies aimed at the detection of circum-stellar 
gas in the AU Mi c disk. By the i r non -detection of far- 
UV H 2 absorption, iRoberge et al. (2005) derived an upper 
limit on the gas column density along the line of sight of 
Nh 2 < 10 19 cm" 2 . The de tection of H.% in fluorescence en- 
abled iFrance et al.l (|2007f ) to derive a column density of 
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N H2 < 2 x 10 17 cm- 2 (T H2 = 800 K 
and 2000 K, respecti vely) . Comparing their H2 value with 
the the CO results of iLiu et al.l (|2004h . they conclude that 
H2 contributes less than about 1/30*' 1 to the total disk 
mass. The H absorption mainly traces interstellar rather 
than circum-stell ar material and h as a column of Ng = 
2.3 x 10 18 cm -2 (jWood et al.ll2005t) . thus corresponding to 
wi thin a factor of five t o the Mg 11 absorption measurements 
of iRedfield fc Linskvl (120021 , N Mg = 1.6 x 10 13 cm" 2 ), as- 
suming solar abundances and Mg II to be the dominant Mg 
species (|Slavin fe Frischll2002D . 

1.3. Disk models 



iKrist et all (|2005f ) used three-dimensional models of the 
scattering cross-section densities throughout the disk to in- 
terpret their HST images. They find that in the inner disk 
region (12 AU < r < 49 AU) the observations can be ex- 
plained by a relatively uniform scattering cross-section den- 
sity (forward scattering particles with an albedo of 0.5) in 
approximate correlati o n with the non-flari n g disk model of 
iMetchev et ail (|2005ft . lAugereau fc BeustJ |2006l) inverted 
the visible and near-IR scattered light profiles to study 
the grain properties and find that the scattered light is 
reflected at grains with sizes mainly between 0.1 /im and 
1 /im. However, they require only about 1/20" 1 of the total 
disk mass to account for the scattered light. 

The dynamics of the grains were included by 
IStrubbe fc Chiangl (|2006l ) into their models and led them 
to explain the observations with a "birth ring" at about 
40 AU, where larger planetesimals of about decimetre size 
are located. By collisions of these planetesimals the sub- 
micron sized grains are produced, which then, depending 
on their size, the radiation pressure, the gas content of 
the disk and the stellar wind, are expelled from t h e sys - 
tem or dragged into AU Mic. IStrubbe fc Chiang! (|2006h 
conclude that AU Mic's debris disk is dominated by de- 
structive grain-grain collisions and that the inner part of 
the disk is largely void of submicron sized grains. These 
small grains mainly populate the outer part of the disk 
and result in the blueish appearance of a debris disk rel- 
ative to the star, since submicron siz ed grains provide th e 
largest fraction of the scattered light (jMetchev et al.ll2005h . 
In contrast to the scattered light, the IR excess is caused 
by larger bodies because the mass of the submicron sized 
grains is too low to a ccount for the obs e rved emission 
([Fitzgerald et al.ll2007t ). lAugereau fc Beustl (|2006l ) already 
noticed the need to increase their disk mass derived from 
the scattered light to reproduce the SED. The void of small 
grains in the inner part of the dis k is co nsistent with the 
polaris ation data of iGraham et all (|2007| ) and the near-IR 
data of lFitzgerald et all (|2007l ). who derived an upper limit 
on the mass of submicron sized grains in the inner zone of 
the disk of 10" 4 M ffi . 

2. The role of X-rays 

The disk models strongly suggest that AU Mic is directly 
shining through its disk, therefore, absorption features re- 
lated to the disk should in principle be present in the ob- 
served spectrum. AU Mic 's strong X-ray emission make 
this system a prime target to search for X-ray absorption 
features from a circum-stellar disk. The dominance of dust 
lets us expect a substantial amount of carbon in AU Mic 's 
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Fig. 1. X-ray transmission curves for different absorbers 
using ba mabs from the PINTofALE-packag e based on the 
data of [Balucinska-Church fc McCammon (|1992[) . Top: 
Absorption by cold gas with solar abundances. Middle 
panel: Gaseous carbon absorption. Bottom: Absorption 
by carbon grains. The individual curves are slightly shifted 
along the x-direction to maintain clarity. The C-K edge is 
located at 43.6 A. 

For the grain absorption, the column density gives the total 
number of carbon atoms along the line of sight. Therefore, 
the number of grains in the line of sight depends on the 
grain-size, i.e., varies between the transmission curves for a 
fixed column density. For comparison, the dashed line in the 
gas absorption panel describes the transmission of carbon 
grains with Nc = 3x 10 18 cm" 2 . 



disk, and thus we examine the influence of carbon on an 
X-ray transmission spectrum. 

2.1. X-ray transmission in the disk of AU Mic 

X-ray transmission spectra contain absorption features di- 
rectly related to individual elements and, therefore, can be 
used to probe the elemental composition of an absorber 
(e.g. iNicastro et ail 120051: IWilliams et"ail l2007t) . In par- 
ticular, X-ray transmission spectroscopy allows an 
assessment of the carbon column density in the disk 
of AU Mic . In Fig. [1] we therefore illustrate the effects 
of X-ray absorption on the transmission curves for differ- 
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ent column densities and the resulting soft X-ray spectrum 
near the carbon K-edge and out to 200 A , i.e., in the 
band pass exclusively accessible to the Chandra LETGS. 
We specifically consider the case of cold gas absorption 
with solar abundances (Fig. [TJ top panel), the case of a 
pure carbon absorber with various column densities (Fig.[TJ 
middle panel), and the case of carbon absorbers locked up 
in various grain sizes (Fig. [TJ bottom panel); clearly, for 
an absorber significantly composed of grains, self-shielding 
within the grains becomes important for the calculation of 
its tr ansmission properties (see Appendix A in lWilms et al.l 
2000). From Fig. [1] (bottom panel) it is clear that graphite 
grains with sizes around 0.1 /jm are ideal to derive column 
densities from the carbon absorption edge, because they are 
so small that self-shielding is unimportant; the derived col- 
umn density represents (almost) the number column den- 
sity of carbon. With increasing grain sizes, the absorption 
features approach that of grey absorbers and grains with 
sizes in excess of 10 /im are virtually completely grey at 
X-ray wavelengths with only a marginal reduction of trans- 
mission. Therefore, the correspondence of absorption edge 
depth and column density holds strictly only for gas and 
small grains, while it breaks down for unknown grain sizes 
(in particular for sizes > 0.3 /mi). 

Fig. [U also illustrates that it is virtually impossible to 
reconstruct the elemental composition of an absorber from 
the reduced transmission outside the absorption edges. The 
shape of the reduced transmission at longer wavelengths 
can be mimicked by choosing appropriate column densi- 
ties for almost every composition, e.g., the differences in 
transmission between a pure oxygen or a pure helium ab- 
sorber in the Chandra LETGS wavelength range longwards 
of 50 A is less than 10 %. Note that only very few narrow 
absorption lines are available in the case of (low ionised) 
disks. In particular the ls-2p transitions of oxygen (23.5 A) 
and carb on (~44.8 A, estimated from the Kg energy) are 
available (jWilms et al.ll200Ct iHenke et al.lll993D . 

In the case of an absorber with solar abundances the 
absorption due to H, He and O reduces strongly the trans- 
mission around the absorption edge of C, however, the 
flux ratio between both sides of the edge depends only on 
the actual carbon column density. For a column density 
of Njj — 10 21 cm -2 assuming cold gas with solar abun- 
dances, the transmission is reduced by a factor of 20, but 
the transmission at the low energy side of the C-K edge 
is about twice that of the high energy side, yet the result- 
ing low flux level makes it virtually impossible to measure 
the depth of the C-K edge. Therefore, using the carbon ab- 
sorption edge requires carbon to be strongly enhanced in 
the absorber so that it contributes significantly to the ab- 
sorption without the overall reduced transmission by other 
elements. 
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Fig. 2. LETGS zeroth order X-ray lightcurve of AU Mic. 
Time represents the time after the start of the observation 
(MJD=54643.506). 



the total mass to small grains, while most of the mass is 
stored in larger bodies with sizes up to ~ 10 cm. In order to 
test these scenarios we calculate the mass column density 
by distributing 0.01 M m uniformly within the inner zone of 
the lKrist et all (l2005h model (10 AU - 50 AU, /i w 2 AU), 
i.e., the zone where the dust causing the infrared excess 
should be located. The mass column density a is given by 



V=-d 



M 



irh(r% — r 
7 x 10~ 4 g/cm 2 



{To 



M 



Trh(r 



(1) 



where M is total disk mass within 50 AU, V is the volume 
occupied by the disk, d is the line of sight and r a and r» are 
the outer and inner radii of the considered region, assuming 
no significant contribution of larger bodies. 

As discussed above, the associated X-ray absorption fea- 
tures depend on the chemical compo sition of the di s k and 
the grain sizes and shapes. Although iRoberge et al.l ([2006) 
found that in the j3 Pic disk the chemical composition might 
deviate from solar by factors of a few with carbon being 
overabundant, we assume for the AU Mic disk that H and 
He are virtually absent and that the mass is stored in the re- 
maining elements with solar abundances. This implies that 
~ 15% of the mass is stored in carbon atoms. The number 
density of carbon atoms in the line of sight is then 



2.2. The mass column density of AU Mic 

The range of possible disk masses and hence column densi- 
ties for AU Mic is relatively narrow. The different methods 
independently point to a total disk mass of 0.01 M m (» 6 x 
10 25 g, cf., Sect. 4.1 in iFrance et ai1l2007t iFitzgerald et al.l 
l2007f ). however, the distribution of this disk mass between 
small grains and larger bodies and their locations are not 
particularly well constrained. The most recent disk models 
including the grain dynamics attribute less than a tenth of 



N C = X, 



c- 



0.15 



m c 



7 x 10~ 4 g/cm 2 
12- 1.66 x 10 - 24 g 



= 5 x 10 18 cm- 2 , 



where Xq is the mass-fraction of carbon (15 %) and mc is 
the mass of a carbon atom. 

From Fig. [T] it is clear that such column densities should 
impose strong feature on the transmission spectrum de- 
tectable with a Chandra LETGS observation. 
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3. Observations, data reduction and immediate 
results 

3.1. Observation and data reduction 

Au Mic was observed on 30 June 2008 with the 
Chandra LETGS (Obs-ID 8894). The total exposure time 
was 50 k sec and the data redu ction was carried out using 
CI AO 4 (jFruscione et alJl2006T ). The photon detector of the 
LETGS is the HRC-S, a micro-channel plate, which does 
not provide sufficient energy resolution to separate the in- 
dividual diffraction orders superposed on the same detector 
area. Therefore, to allo w an analysis o f the data with stan- 
dard tools like XSPEC C& rnau d 1996). we constructed new 
response matrices including up to ninth diffraction order 
contributions following the instructions given in the CIAO 
thread^]] (see Appendix [A")) . Line fluxes were obtained using 
CORA which accounts for the Poisson character of the data 
(iNess fe Wichmannl[2002h . 

We experimented with the "standard" filters (light, 
medium and heavy) and with the new Gain Map and 
Pulse-Height filteiQ Around the carbon edge, the differ- 
ence in the background fraction is only a few percent be- 
tween the light and medium filter while significant for the 
new filter; we list the corresponding values in Tabled) The 
figures and numbers given in the text pertain to the stan- 
dard light filter, unless otherwise noted. 



3.2. Global plasma properties 

In Fig. [5] we show the zeroth order X-ray count rate of 
AU Mic vs. time; obviously, the light curve is more or less 
constant during the first 40 ksec of the observations, after- 
wards there is a small flare-like increase. We do not treat 
this increase separately but analyse the observation in to- 
tal. In Fig. [3] we show the recorded LETGS spectrum up to 
120 A; sub -regions will be shown individually in the next 
sections. AU Mic 's X-ray spectrum is typical of an active 
star, showing the strong Ne emission features and a strong 
OVIII Ly a line; the carbon Ly a line is also quite strong. 
The OVII triplet between 21.6 and 22.1 A as well the Fe 
XVII feature at 15.03 A are relatively weak as well as all 
features attributable to N; we do point out the emission 
lines between 90 and 120 A attributable to highly ionised 
iron. We used XSPEC to fit the X-ray spectrum using a 
combination of th ree APED models (variable abundances, 
ISmith et al.ll200ll) with one absorption component. The in- 
creasing background at longer wavelengths make these bins 
less useful for the spectral analysis. We therefore restrict the 
wavelength range to values between 5 A and 35 A. Using 
the full wavelength range does not influence the two low 
temperature components; only the best fit temperature of 
the high temperature component doubles. This is mainly 
driven by a single, somewhat strangely shaped Fe XX line 
at 132.85 A. Furthermore, the abundances of Mg, Al, Ca 
and Ni (low FIP) were coupled to that of Fe in order to 
decrease the number of free parameters. The thus obtained 
fit results (listed in Tab. [D show the inverse FIP effect 
usually found in M-dwarfs piobrade fe Schmittl[2005h . The 
temperature structure and the abundances compare well 
with a fit performed with the XMM- Newton RGS-data of 
AU Mic (jobs-id 01114 20101'). 



Table 1. Coronal proper ties of AU Mic 
given relative to that of lAsplund et al.l (|2009t) and 90 
errors are given. 



Abu ndances are 

% 



Property 



Value 



N H 


2 x 10 is cm" 2 


feTi 


0.29 ± 0.02 keV 


EMi 


7.5+J g x 10 51 cm 5 


kT 2 


0.67±°'°3 keV 


EM 2 


11.0+|g x 10 51 cm 5 


kT 3 


1.49lg;S keV 




5.6+i'l x 10 51 cm 5 



Element 


Abundance 


C 


0.9 ±0.3 


N 


0.8 ±0.3 


O 


0.5 ±0.2 


Ne 


-i q+0.3 
1 " a -0.1 


Si 


0.3 ±0.1 


s 


1+ ' 2 


Fe 


n IQ+006 
u - iy -0.04 



1 http : //exc .harvard. edu/ciao/threads/hrcsletg_orders/ 

2 http : / / exc .harvard. edu/ contrib/letg/GainFilter/ 



4. Absorption at the carbon edge 

The carbon K-edge at 43.6 A (284 eV, i.e., the en- 
ergy needed to expell a K-shell electron from an isolated 
carbon atom) is a promising absorption feature in the 
Chandra LETGS wavelength band. In this region, the line 
emission is relatively small compared to continuum emis- 
sion. The oxygen edge, located at 23.1 A (536 eV), is also 
a good candidate, but with an enhanced line-to-continuum 
ratio. This increases the errors since the prediction of the 
line emission is not possible with the required precision. 
Unfortunately, the UV/Ion shield of the HRC-S also con- 
tains a large amount of carbon (Nc ~ 10 18 cm -2 ) leading 
to an instrumental feature very similar to that of expected 
interstellar/circum-stellar carbon absorption. Still, the high 
continuum fraction of the first order emission around the 
C-K edge makes this region more promising than the O-K 
edge. 

4.1. The model 

As is obvious from Fig. [T] the height of the C edge does 
depend sensitively on the carbon column density (and on 
the grain sizes). Therefore our basic idea to measure the 
height of that edge is as follows: We divide the wavelength 
range around the C edge into a low (44 - 48 A) and a high 
energy band (30 - 43 A) and sum up the counts in these 
two bands to increase the signal. By comparing these two 
count numbers with models containing specific amounts of 
carbon absorption the best fit carbon column density can 
be found. 

We do not extend the wavelength range to longer wave- 
lengths since these longer wavelength bins are close to 
the detector gap at 50 A of the negative diffraction or- 
der, where the data is less reliable. We also exclude the 
immediate region around the edge and around the ls- 
2p carbon absorption line (^44.8 A) because of the rela- 
tively little known fine-structure of the carbon absorption 
edge. Exclusion of a larger wavelength region around the 
edge does not change the results significantly (see Fig. [7]). 
Increasing the amount of carbon in the line of sight, changes 
the flux-ratio between the low and high band; since the the 
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Fig. 3. The measured spectrum in the spectral range relevant for the analysis. Shown is the sum of the positive and the 
negative diffraction order. 
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Fig. 4. Model and data around the C-K edge (sum of both 
diffraction orders). 



change in transmission in the low energy band is compara- 
bly small, we normalise our model to match the low energy 
band. Therefore we carry out only a differential analysis. 



Table 2. Contributions at the C-K edge (summed contin- 
uum and line emission) . Listed are the total counts of the 
complete model (sum of positive and negative diffraction or- 
der). Top is the result from the standard light filtering and 
in the bottom the corresponding value for the new Gain 
Map and Pulse-Height filtering are shown. 



Contribution 




Wavelen 


gth range 






30 


- 43 A 


44 - 


48 A 


Light filter: 

Total counts 


6495 




2593 




First order emission 


1546 


23.8 % 


1014 


39.1 % 


Higher order contribution 


569 


8.8 % 


172 


6.6 % 


Background 


4251 


65.4 % 


1287 


49.6 % 


Gain Map and Pulse-Height filter: 






Total counts 


4582 




1939 




First order emission 


1467 


32.0 % 


949 


48.9 % 


Higher order contribution 


565 


12.3 % 


173 


8.9 % 


Background 


2376 


51.9 % 


705 


36.4 % 



For our modelling of the recorded spectrum above and 
below the carbon edge we consider the following four com- 
ponents as contributors to the counts around the C-K edge. 

1. The continuum emission. 
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Fig. 5. The different contributions to the measured counts 
around the C-K edge. Top: First order emission, middle: 
Higher order contribution, bottom: Background. 

2. The superimposed line emission. 

3. The higher diffraction orders which cannot be filtered 
out due to the low intrinsic energy resolution of the 
photon detector (HRC-S). 

4. The enhanced background of the LETGS due to a wiring 
error of the micro-channel plate detector. 

In the following we describe our modelling of the individ- 
ual components and discuss their accuracy; the given errors 
pertain to the higher energy side of the edge while the num- 
bers in brackets correspond to the lower energy side. Since 
the models are normalised to match the low energy side of 
the edge, an error at this side will also influence the num- 
ber of predicted counts at the high energy side of the edge 
(by roughly the same amount). Fig. H] shows the complete 
model and the individual model contributions are shown in 
Fig. O 

The continuum is the simplest component since its 
shape around the C-K edge does not change noticeably 
for reasonable changes in the plasma properties. The main 
spectral components are known from our global spectral 
modelling, but changing, for example, the temperature by 
about 50 %, changes the number of counts left of the carbon 
edge by less than 2 % (50 cts), while preserving the counts 
at the right side of the edge by choosing an appropriate 
normalisation. 

The treatment of emission lines is less straightforward. 
We estimate that contribution of the first diffraction order 
emission lines is not predictable with a better than 10 % 
accuracy. Therefore, a region of 0.2 A around each strong 
line (based on our XSPEC / APED emission model) was ex- 
cluded from the two energy bands thus removing ~ 98% of 
the line flux of these strong lines. Only 65 % and 73 % of 
the available bins are therefore used in the 30 - 43 A and 
44 - 48 A range, respectively. As a consequence, the ratio 
between line and continuum emission is below 10 % in the 
high and low band. In addition, the remaining counts from 
unresolved lines are predicted from the emission model and 
are added to the model as a correction; these numbers turn 
out to be 50 and 35, respectively. Taking this numbers as 
a conservative estimate of the error, the impact of these 
lines on the derived carbon column density is quite small, 
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Fig. 6. Data and model for the sum of the positive and 
the negative diffraction order. Note that the model is not 
based on a physical model (except for the continuum and 
the absorption). 



therefore, also the influence of " unknown" lines ought to be 
small. 

The third diffraction order of the LETGS is the 
strongest higher order contribution to the region around the 
C-K edge, since the second diffraction order is suppressed 
by the grating design. Unfortunately, AU Mic being an ac- 
tive star, the Ne IX and Ne X lines are quite strong, and 
their third order components are easily discernible in our 
LETGS spectrum (cf. middle panel in Fig. [5]). However, first 
and third order are always measured simultaneously, thus, 
by fitting the stronger (identifiable) lines independently in 
the wavelength range shorter than the C-K edge provides 
the desired description of the higher order contamination. 
The selection of the lines is based on the APED model (see 
Sect. O, the Chandra HETGS spectrum of AU Mic and 
on the lines visible in the Capella spectrum (Capella is 
a calibration target of the LETGS, thus providing about 
400 ks of data and has an approximately comparable tem- 
perature structure as AU Mic). Note that this method is un- 
affected by the uncertainties of the spectral model. Figure [5] 
shows the model of the wavelength range which provides 
the strongest contamination at the carbon edge; 97 % of the 
photons in this range are included in the model. Taking the 
remaining 3 % as a measure for the accuracy of the higher 
order model, these 3 % correspond to only 12 (4) photons 
not contained in the higher order model between 30 A and 
48 A. We can therefore neglect this effect on the derived 
carbon column density compared to the other factors like 
the statistical noise or a potential error in the higher order 
diffraction efficiencies, which would for a 10 % error0 result 
in an error on the 30 count level (6 cts). 

The last component to consider is the background. 
Fortunately, the background is relatively uniform in the 
wavelength range between 30 A and 48 A. Since the bins 
in the "line-free" regions are all summed up, the impact 
of the short-scale variation should be further reduced. The 
largest problem of the background in this analysis is the 
enlarged count number leading to an increased statistical 
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uncertainty which amounts to about 50 (30) counts. Since 
the models are tuned to match the longer wavelength side 
of the edge, the combined statistical error is ~ 75 counts 
at the high energy side of the edge. 

4.2. The influence of carbon absorption 

Varying the absorbing carbon column density changes 
the transmission curve in the whole wavelength range 
(see Fig. [TJ. Therefore, it is necessary to re-adjust the 
continuum-level and the amplitudes of the emission lines 
for the higher order contribution to construct a new best 
fit model. The influence of the carbon absorption on the 
low-energy side of the edge is relatively small (cf. Fig. [1]), 
therefore, we tuned the models so that this region matches 
the data by varying the continuum normalisation of the 
model. The normalisations of the models with zero absorp- 
tion and N c = 2 x 10 18 cm" 2 differ by about 20 %. Due 
to the strong line emission (there are not many wavelength 
regions without line emission at shorter wavelengths) and 
the refitting of these lines, this does not noticeably change 
the fit-quality at the shorter wavelength. Having thus fixed 
our model to account for the number of counts on the low- 
energy side of the carbon edge, we can compute the model 
counts on the high energy side of the carbon edge as a func- 
tion of the assumed carbon column density and hence, by 
comparison with the observation, the missing model counts, 
since increasing the carbon column density leads to a more 
and more relatively reduced model count number below the 
edge. The dependence of this number of missing counts in 
the model on the assumed carbon absorption is graphically 
displayed in Fig. [JJ Depending on the accepted accuracy 
of our model we can read off the maximum number of car- 
bon atoms along the line of sight assuming gaseous ab- 
sorption (i.e. no self-shielding). If grains cause the absorp- 
tion, a larger number of carbon atoms in the line of sight 
is required to produce the same number of missing counts. 
Explaining gaseous absorption with Nc — 10 18 cm -2 re- 
quires a 1.3 times higher number of carbon atoms in the 
line of sight for 0.1 fjxa sized grains, a 1.9 times higher 
number for 0.3 ^m and an almost 5 times higher number 
for 1.0 /jm grains. 

Collecting the errors attached to the individual model 
components, we estimate an overall error of approximately 
200 (80) counts. We cannot simply add the error of both 
sides of the edge, since it is rather unlikely that the same 
error-source results in an erroneously increased number of 
counts at one side of the edge while, at the same time, 
underestimating the counts at the other side of the edge. 
Therefore, we conclude that a robust estimate of the error 
is about 200 counts or Nc < 10 18 cm~ 2 (for pure carbon 
gas). 

The result shown in Fig. [JJ suggests only a small amount 
of carbon-atoms along the line of sight (Nc ~ 10 17 cm~ 2 ), 
which might be supported by a slight deficit in counts 
around the carbon ls-2p absorption line. However, this is 
probably only a statistical effect since this effect is reduced 
by using the stronger Gain Map and Pulse-Height back- 
ground filtering (see dashed line in Fig. [JJ). 

5. Constraining the total absorption 

The method described above offers the ability to directly 
constrain the column density of carbon in the line of sight. 



400 F 



o 




-100 E 



0.0 0.5 1.0 1.5 2.0 

Carbon column density (10 18 cm -2 ) 

Fig. 7. Number of missing (positive count number) or over- 
predicted photons (negative count number) as a function 
of the assumed carbon absorption column density in the 
model. The flattening of the curves indicates that the flux 
will, eventually, be completely absorbed in the considered 
wavelengths range for even higher column densities. The 
dashed line is the average of the positive and negative or- 
der using the new Gain Map and Pulse-Height filtering. 
The dashed-dotted line is the result using the wavelength 
ranges 30-42 A and 45-48 A. The statistical errors reflects 
only first order errors. For a discussion of the errors see 
Sect: gj 

It is, on the other hand, not useful to set tight contrains on 
the "total" absorbing column density since only a narrow 
wavelength range is inspected. Emission lines located at 
largely different wavelengths better utilise the large wave- 
length coverage of the LETGS. Usage of a single element 
and ionisation stage reduces the dependence on the recon- 
struction of the abundances and the temperature structure 
of the corona. 

A good candidate for such a study is Fe XVIII with 
lines at 14.21 A (blend of 14.2060 and 14.2085 A), 
16.08 A (blend of 16.0760 and 16.0913 A), 93.92 A and 
103.94 A. Assuming that all these lines are produced in 
the same environment, the relative fluxes of these ions de- 
pend only weakly on the temperature structure; the de- 
pendence on the density is virtually negligible for densi- 
ties around n K ~ 10 10 cm" 3 as usually found in stellar 
coronae (e.g. iNess et al.ll2002l ). Unfortunately the stronger 
Fe XVIII lines at shorter wavelength are blended. While 
the 14.21 A lines do posses only small contributions 
from the 14.17 A Fe XXI line and the Fe XVIII doublet 
at 14.26 A , the 16.08 A is blended by the strong 
O VIII doublet at 16.00 A and a Fe XIX line at 16.11 A, 
which amounts to about a fourth of the Fe XVIII emis- 
sion for the temperatures in question. The derived fluxes of 
the short wavelength lines are comparable to the available 
HETGS data (Obs-ID 17). 

The Chianti package (jDere et al.l Il997t lLandi et al.l 

I2006h was used to predict relative line fluxes. To investi- 
gate the temperature dependence of the relative fluxes, a 
Gaussian fit of the emission measure distribution (EMD) 
of Appendix [C] was performed. The line ratios were then 
calculated for Gaussian EMDs with different peak temper- 
atures and widths, e.g. the peak temperature was changed 
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Fe XVIII 93.92 A 



Fe XVIII 103.94 A 



Fe XVIII 16.08 A 



6.8 7.0 

log (T/K) 



Fig. 8. Ratio of the line-fluxes to the Fe XVIII 14.21 A line 
using Gaussian-like emission measure distributions with a 
width of log(T/K) = 0.5. 



up to Alog(T) = 0.4 (see Fig. [5]). The measured fluxes 
(sum of both LETGS orders) of the above lines together 
with predicted ratios (the 14.21 A line normalised to 1.0) 
are listed in Tab. [3] 

From the difference between predicted and measured 
counts the required hydrogen column density (assuming 
solar abundances) can be calculated. The resulting ranges 
are listed in Tab. [3J they include the uncertainty in the 
reconstruction of the temperature and the statistical er- 
ror. The influence of the adopted temperature structure 
is small compared to the statistical errors. Comparison 
of the line fluxes from the standard filtering and that of 
the new Gain Map and Pulse-Height filtering shows that 
both values agree well. However, the upper limits from 
these fluxes are higher than that from the light filtered 
data since their line fluxes at longer wavelength are lower 
(N H < 11 x 10 18 cm- 2 ). 

Summing the line fluxes of the short wavelengths lines 
and combining the upper limits of the two long wavelength 
lines improves the upper limit to Njj < 4 x 10 18 cm~ 2 (Ict). 
Instead of using one single ion, we checked the ratio 
of the two strongest Fe lines at short and long wave- 
lengths (Fe XVII 15.02 A and Fe XX 132.85 A). Their ra- 
tio depends more strongly on the temperature structure 
with predicted ratios (132.85 A/15.02 A) between 1.0 
(log T=6.5)and 2.0 (log T=7.3) for a shallow EMD (AT = 
0.5). However, their resulting upper limit is only slightly 
higher (N H < 6 x 10 18 cm" 2 ) than that from the Fe XVIII 
lines. 

The X-ray absorption at wavelengths of around 100 A is 
caused mainly by He atoms. The upper limit on He from 
the observed line fluxes is Nn e < 6 x 10 18 cm -2 assuming 
that the absorption is caused exclusively by He and Nh < 
14 x 10 18 cm~ 2 for a pure Hydrogen absorber. Similarly, 
upper limits on the abundance of other absorbing elements 
can be derived under the assumption that they are the only 
absorbing species. 

The analysi s of the absorption towards Capella by 
IGu et all (|200fift using a comparable method has an error 
of about a factor of five lower which is consistent with the 
larger data base available for Capella which is a calibration 
target of the Chandra LETGS. 



6. Discussion 

Neither the edge based method nor the line based method 
are able to find significant absorption along the line of 
sight towards AU Mic. For both methods, the statistical 
error overwhelms the error caused by the incomplete mod- 
els. This is clear by inspection of Tab. |3]for the line based 
method but is also true for the edge analysis. The detectable 
carbon column density from an analysis of the edge height 
is 5 ... 50 x 10 17 cm -2 for artificial stars with X-ray fluxes 
within a range of five around that of AU Mic (50 ks LETGS 
exposure) . This value decreases only slowly for even higher 
X-ray fluxes and "saturates" at Nc ~ 3 x 10 17 cm~ 2 due 
to the potential error in the higher order contribution (see 
Appendix [B]). The edge-method provides a direct upper 
limit on the carbon content of the disk without any as- 
sumption on the absorber. Its limit is about that provided 
by the line based method assuming that only carbon is in 
disk, i.e., calculating the maximum carbon column density 
allowing the detection of the line-flux from the long wave- 
lengths Fe XVIII lines. 

Two measurements of H (atomic and molecular) re- 
strict the total hydrogen column density in the line of 
sight towards AU M ic to Nh <4x 10 19 cm~ 2 (mainly H 2 , 
iRoberge et al.ll2005h . From this value, on ly 2 x 10 18 cm~ 2 
are atomic hydrogen ([Wood et al.l I2005D . Note that the 
conversion of mol ecular hydrogen t o Nh in stellar disks 
is uncertain (e.g., iKamp et al.l l2007h . while molecular hy- 
drogen is orders of magnitude less abundant than atomic 
hydrogen in the n earby interstellar medium (d < 200 pc; 
iLehner et al.l 2003T) . Th e locat ion of the fluorescent H2 de- 
tected by iFrance et al.l (|2007f ) is not necessarily along the 
line of sight and leaves space for an additional H2 not con- 
tributing to the fluorescent H2 emission. 

X-rays are sensitive to both, atomic and molecular, hy- 
drogen. H2 absorption i s about 2.8 times stronger than 
that of a single H atom (jWilms et al.ll20 00). Furthermore, 
the X-ray absorption is insensitive to the excitation state 
of i?2- Therefore, the ha upper limit from the UV mea- 
surement (Nh 2 < 17 x 10 18 cm -2 ) corresponds to about 
Nh ;$ 50 x 10 18 cm~ 2 in soft X-rays; the same order as 
the upper limit derived from the line fluxes assuming a 
pure hydrogen absorber. The X-ray upper limit on H also 
includes H2 in the line of sight located in the inner and 
outer part of the disk which might have different excita- 
tion stages and ther efore complements the upper limit of 
IRoberge et al.l (|2005l ) . Furthermore, the detected line fluxes 
at A > 100 A can be translated to upper limits on other 
elements restricting the gas / small grain content of the 
disk, e.g., that of carbon in gaseous form and small grains 
(s < 0.3 fim, see Tab.©. 



7. Summary 

We analysed the impact of absorption caused by AU Mic's 
debris disk on its observed X-ray spectrum resulting in 
three upper limits on the column densities for hydrogen, 
helium and carbon: 



1 1- 2 (p Ure jj absorption) 
2 (p Ure jj e absorption) 



- N He < 1.2 x 10 

— Nc < 10 18 cm~ 2 (pure C absorption) 

Assuming an absorber with solar abundances, the upper 
limit is Nh < 10 19 cm~ 2 , which is about a factor of five 
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Table 3. Strong Fe XVIII lines and measured fluxes. The values for the new Gain Map and Pulse-Height filtering 
are also given. Their corresponding ratios are given in brackets in the ratio column while the values in the columns for 
the corresponding absorbing column densities are based on the standard filter and the value in brackets give the upper 
limit including the statistical error. 



Wavelength 


Measured flux Measured flux^ 


Ratio to 


Predicted 


Corresponding 


Corresponding 




(10" 5 ph/cm7s) 


14.21 A 




N H (10 1S cm" 2 ) 


N c (10 17 cm" 2 ) 


14.21 


6.3 ± 1.4 6.3 ± 1.4 










16.08 


5.7 ± 1.3 6.3 ± 1.0 


0.9 ± 0.3 (1.0± 0.3) 


0.6^ 






93.92 


28.9 ± 3.7 25.8 ± 3.3 


4.6 ± 1.1 (4.1 ± 1.1) 


3.9 ±~0.3 


0-7 


0-8 


103.94 


11.1 ± 3.4 9.7 ±2.7 


1.8 ± 0.7 (1.5 ± 0.5) 


1.4 ±0.1 


0-8 


0-10 



a New Gain Map and Pulse-Height filter 
b plus 0.2 blend 



higher than the pure interstellar value. The statistical er- 
ror caused by counting statistics and the high background 
prevents to set lower limits, while the presented methods 
are in principle more sensitive. 

Both upper limits are consistent with the idea that the 
AU Mic disk is optically thin in the radial direction as pro- 
posed by recent analyses from optical and infrared measure- 
ments, and in line with current disk models, which assume 
that the inner part of the disk is almost void of small grains. 
The debris of the collisions of planetesimals in the "birth- 
ring" populate mainly the outer part of the disk, where 
their density is so low that they escape a detection in this 
X-ray observation. In the inner part of the disk larger grains 
hold the mass as predicted by birth-ring scenarios. 

iRoberee et al.l (|2005l ) state that they find weak signs 
of H2 absorption on the order of 10 18 cm" 2 in their UV 
data. To reach this level with the Chandra LETGS setup, 
longer integration times are required. For Capella, a cali- 
bration target of the LETGS, the N H - 10 18 cm" 2 range 
is accessable. However, it seems unlikely to reach down 
to Nh ~ 10 17 cm" 2 with the present instrumentation 
as would be required to safely distinguish a hypothetical 
circum-stellar H2 component in the Nh ~ 10 18 cm" 2 range 
from the interstellar atomic H absorption. 

Appendix A: Higher order response matrices for 
XSPEC 

The response matrices contain the probabilities describing 
how the detector will respond on a photon with a given 
energy. Including the higher diffraction orders therefore re- 
quires to add additional probabilities at two, or more times 
the rest wavelength. Since the resolution grows with the or- 
der number, the energy grid needs to be refined accordingly. 
The FITS-specification for response files (CAL/GEN/92- 
002) already includes the concept of wavelength groups 
which are ideal to reduce the size of the final response ma- 
trix. 



Appendix B: Carbon detectability 

To estimate the detectability of carbon absorption from the 
edge height, we need to quatify the different error contri- 
butions as a function of the source count-rate. Concerning 
the error due to unknown unresolved first order emission 
lines and uncertainties in the temperature structure, we re- 
gard a correlation with the statistical error as realistic, since 



better data quality enables a more precise temperature re- 
construction and, in turn, a better model for the emission 
lines around the carbon edge. However, the detailed de- 
pendence on the source count-rate might be more compli- 
cated. The higher order contribution scales with the source 
flux; therefore, its error scales also linearly with the first or- 
der flux since calibration uncertainties dominate the error. 
Figure IB. II shows the minimum detectable carbon column 
density for different source fluxes assuming a pure carbon 
disk. It shows that column densities below 3 x 10 17 cm" 2 
can only be achieved for star an order of magnitude X- 
ray brighter than AU Mic and with a hypothetical 500 ks 
observation. However, the existence of other elements in 
the disk increases the detectable column density due to the 
absorption of the first order flux around the carbon edge. 
Furthermore, without knowing the composition of the disk, 
the extra absorption changes the continuum slope around 
the carbon edge and therefore increased the uncertainty in 
the carbon column density. 



10 1 



E 







AU Mic 










500 


S"^^ — 
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Fig. B.l. The detectable carbon column density by mea- 
suring the jump height at the carbon edge as a function of 
the X-ray flux (0.1 keV - 2.0 keV). Here, a pure carbon disk 
is assumed. The two lines represent a 50 ks and a 500 ks ex- 
posure with the LETGS. The same temperature structure 
as AU Mic is assumed which determines the higher order 
contribution around the carbon edge. The vertical dotted 
line indicates the X-ray flux of AU Mic. 
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Appendix C: The emission measure distribution 

The large wavelength coverage of the LETGS offers the abil- 
ity to quantify the emission measure distribution (EMD) 
in the corona. Using the lines listed in Tab. IC.f 1 we recon- 
structed the EMD using the PintOfAle (Kashva p fe Drake! 
2000) package. Our result compares well the result obtained 
usi ng the line fluxes measu red in the UV wavelength range 
by iDel Zanna et al. | (120021) . We als o experimented wit h the 
Chianti package (jDere et al.l Il997t lLandi et al.l [2006) and 
found that the results are strongly influenced by numerical 
problems due to the shallow temperature coverage of the 
observed lines. However, both methods give rather similar 
results. 



Table C.l. Strong lines and measured fluxes. 




Ions 


Wavelength 


Flux (10- 5 ph/cm7s) 


Si XIII 


6.65 


4.5 ± 0.7 


Ne X 


10.24 


6.3 ± 1.0 


Ne IX 


11.54 


6.3 ± 1.0 


Ne X 


12.13 


44.1 ± 2.1 


Ne IX 


13.45 


32.8 ± 1.9 


Ne IX 


13.54 


11.2 ± 1.3 


Ne IX 


13.70 


22.9 ± 1.6 


Fe XVIII 


14.21 


6.3 ± 1.4 


Fe XVII 


15.01 


18.0 ± 1.4 


Fe XIX/O VIII 


15.18 


6.4 ± 1.1 


Fe XVII/Fe XX 


15.26 


8.4 ± 1.1 


O VIII 


16.01 


17.6 ± 1.5 


Fe XVIII 


16.08 


5.7 ± 1.3 


Fe XVII 


16.78 


12.1 ± 1.2 


Fe XVII 


17.10 


30.2 ± 1.9 


O VIII 


18.97 


110.3 ± 3.3 


O VII 


21.61 


29.2 ± 2.3 


O VII 


21.79 


10.0 ± 1.6 


O VII 


22.10 


19.4 ± 2.0 


N VII 


24.79 


16.6 ± 2.0 


C VI 


33.73 


33.9 ± 2.9 


Fe XVIII 


93.89 


28.9 ± 3.7 


Fe XVIII 


103.92 


11.1 ± 3.4 


Fe XIX 


108.33 


26.2 ± 4.3 


Fe XX 


118.66 


14.1 ± 4.2 


Fe XX 


121.85 


17.0 ± 5.0 


Fe XXI 


128.75 


<16.2 


Fe XX 


132.85 


37.7 ± 8.6 


Fe XXII 


135.78 


27.6 ± 7.9 



Fig. C.l. The reconstructed emission measure distribution. 
The dotted line indicates the Gaussian fit to the EMD. 

It is clear, that most of the emission is produced in the 
temperature range 6.4 < log(T/K) < 7.2. We also show in 
Fig. IC.ll the location of the three temperature components 
from the XSPEC fit with APEC models (variable abun- 
dances). For an estimate of the error in the relative line- 
fluxes in Sect. [S] we fitted the emission measure distribution 
with a single Gaussian. We can then change the centroid 
and the width of the Gaussian to analyse the impact of an 
error in the reconstructed temperature distribution on the 
relative linefluxes. 
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